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Abstract 

Multiple 3d orbitals contribute to the low energy electronic structure in the newly discovered 
iron-based superconductors. The specific orbital characters of the multi-bands are crucial for 
constructing microscopic models. We here report strong polarization dependence of the angle- 
resolved photoemission spectra of BaFei.85Coo.i5As2. Our analysis identifies the orbital characters 
of the low lying electronic structure to be quite different from the existing calculations. 

PACS numbers: 74.25.Jb,74.70.-b,79.60.-i,71.20.-b 
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The discovery of high temperature superconductivity in iron pnictides has ignited inten- 
sive studies on these materials [l-3|. In contrast to the cuprates, where only the dx^i-y^ 
orbital out of the five Cu 3(i-orbitals is important in the low lying electronic structure, the 
iron-based superconductors exhibit complex multi-band nature. Density functional theory 
(DFT) calculations show that its electronic structure near the Fermi energy (Ep) is mainly 
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consisted of all the five Fe 3d orbitals 14, l5ll. Further analysis shows that it is necessary 
to involve at least four bands in the Hamiltonian to model the band structure from DFT 
calculations g], 7|. However, there have been inconsistencies over the orbital identities of the 



bands amongst theories j4-7|, and a few scattered experiments 8l-ll l| . The orbital characters 
are directly related to physical properties in both spin density wave (SDW) states [12j and 
superconducting states |13l4l5l|. It is thus essential to experimentally identify the orbital 
characters of the bands near Ep, which would facilitate the construction of correct models 
for the iron-based superconductors. 

In this paper, we report a systematic study of the orbital nature of the electronic structure 
of BaFei.85Coo.i5As2 near Ep by exploiting the polarization-sensitivity of the orbitals in 
angle-resolved photoemission spectroscopy (ARPES). We have observed strong polarization 
dependence for all the bands near Ep, and identify their possible orbital characters by 
analyzing the intensity asymmetry detected on opposite sides of F. The observed polarization 
dependence can not be understood in the current available five-band models jo^T and also 
differs from that in the SDW states of parent compounds reported recently 9 
result gives strong constraint on the theory of iron-based superconductors and suggests 
that the doping may strongly influence the band structure due to strong electron-electron 
correlations. 

BaFei.85Coo.i5As2 single crystals were synthesized by self-fiux method [l^ with a super- 
conducting Tc of 25K, and no SDW or structure transitions. As a result, the low temperature 
electronic structure is free from the complications of the band splitting and folding 
Data were taken with 40 eV photons at the Beamline 1 of Hiroshima synchrotron radiation 
center (HSRC), or with 100 eV photons at the SIS beamline of Swiss Light Source (SLS). 
By rotating the light polarization at SLS or the entire photoemission spectroscope around 
the incoming photon beam at HSRC, two different polarization geometries were achieved as 
shown in Fig. [T](a). The incident beam and the sample surface normal define a mirror plane. 
For the n (or a) experimental geometry, the electric field direction of the incident photons 
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FIG. 1: (color online) (a) Experimental setup for polarization-dependent ARPES. (b) Photoemis- 
sion intensity map of BaFei.85Coo.i5As2 at Ep in the Brillouin zone (dashed line), kj; and ky are set 
with the unfolded Brillouin zone (solid line). The data were taken in the it geometry with 100 eV 
photons. The measured Fermi surfaces are shown in red lines, among which the Fermi pockets 
around M induced by zone folding are plotted in dashed lines. 

is in (or out of) the mirror plane. All the data were taken with Scienta electron analyzers 
at lOK, the overall energy resolution is 15meV, and angular resolution is 0.3 degree. The 
sample was cleaved in situ, and measured under ultra-high-vacuum of 1 x lO^^^torr. 

Photoemission intensity map at Ep for BaFci 85C00.15AS2 is shown in Fig. [T](b). There are 
two hole pockets around F and several electron pockets around M. Because of the appearance 
of Arsenic ions, there are two iron ions per unit cell. The unfolded Brillouin zone for one iron 
ion per unit cell is shown by solid lines in Fig. [l](b). The electron pockets near M plotted in 
dashed lines are caused by zone folding. We assign F — M as the kx direction. 

Let us first consider the electronic structure around F in the F-M direction. Photoemission 
data taken in a and vr geometries are compared in Fig.[2l Two sets of bands could be observed 
in the vr geometry (assigned as a and /3 band respectively), while only the a band is obviously 
observed in the a geometry (there are some very weak traces of /3 band near Ep). Note that 
the dispersion of the a band changes slightly in two experimental geometries [Figs. [2](c) and 
Et^g)], possibly due to the influence of the (3 band in the tt geometry. 

The matrix element of the photoemission process can be described by 

\M%\ ^ \{ct>)\e ■ r\<pf)\' 

, where e is the unit vector of the electric field of the light For high kinetic-energy 

photoelectrons, the final-state wavefunction 0^ can be approximated by a plane-wave state 
e^^'^ with k in the mirror plane for the analyzer slit setup in Fig. [T]^a). As a result, it is even 
with respect to the mirror plane. For the photoemission process to be symmetry-allowed, 
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FIG. 2: (color online) The polarization dependent APRES data around T along the F-M direction, 
(a) The momentum cut for the presented data is marked with a thick solid line in the Brillouin 
zone. The symmetries of the five 3d-orbitals are illustrated, (b) Photoemission intensity taken 
with IT geometry, (c) The second derivative with respect to energy of the data in panel b, whose 
minimal locus gives the band structure visually, (d) The energy distribution curves (EDO's) for 
the data in panel b. (e) The momentum distribution curves (MDC's) and the fitting results for 
the data in panel b. (f, g, h and i) are the same as in panel b, c, d and e respectively, but taken 
with a geometry. 

{i ■ r)\(f)f) must be even. In turn, it gives that the initial state wavefunction must be 
even (or odd) with respect to the mirror plane in the vr (or a) experimental geometry. 

With respect to the xz-plane, the dyz and dxy orbitals are odd, and the dxz, dx^^y2 and d^^ 
orbitals are even, as illustrated in Fig. [2]^a). Therefore, the /3 band that almost disappears 
in a geometry is dominated by the d^z, d^^-y-i or dz^ orbitals. On the other hand, the a 
band could be observed in both geometries. Based on the MDC's in Figs. [2](e) and[2|^i), the 
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FIG. 3: (color online) The polarization dependent APRES data around M along the F-M direction, 
as marked with a thick solid line in the Brillouin zone in (a), (b) The sketch of the band structure 
in panel d and g. (c) Photoemission intensity and (d) its second derivative with respect to energy 
taken with vr geometry along the cut shown in panel a. (e) The MDC's for the data in panel c. (f, 
g and h) are the same as in panel c, d and e, but taken with a geometry. Data were taken with 
100 eV photons. 



a band has even stronger intensity than the (5 band in vr geometry. Therefore, the a band 
contains both even and odd orbitals significantly. 

Further along the F — M cut, the states around M also show strong polarization depen- 
dence [Fig. 13. Band calculations predicted two bands around M that cross Ep as sketched 
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in Fig. [SJl^b) [4, |5|, |7[, which have been observed before [9|, [ICj]. The 7' band shows up in the 
vr geometry [Figs. [3](c-e)], while the 7 and a' bands show up in the a geometry [Figs. |3](f-h)]. 
We can thus conclude that near Ep, 7' is primarily consisted of dxz, c^x^-^a, or d^-i orbitals, 
while 7 and a' are made of d^y or dyz orbitals [Fig. El^g)]. It has been suggested that both 7 
and 7' bands are of even symmetry in the SDW state of SrFe2As2 [9] . However, note that our 
data are consistent with the normal state data of BaFe2 As2 in Ref . [lOj , thus, the discrepancy 
suggests that the SDW or doping might influence the electronic structure significantly. We 
leave this to further investigations in the future. 

The observed orbital characters are hard to be understood within the current band cal- 
culation results 0, ~\. Considering the five-band models constructed by fitting the LDA 
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FIG. 4: (a), (b) and (c) are the simulation of the matrix element of dxz, and dx2_y2 in vr 
geometry respectively, (d) and (e) are the simulation of the matrix element of dyz and dxy in a 
geometry respectively. The arrows in panels a and b indicate the momentum position where the 
experimental intensity asymmetry can be achieved. 



results in the unfolded Brillouin zone, the two hole pockets around F mainly contain dxz 
and dyz orbitals. More explicitly, along F — M, states at the inner (a) and outer (/3) pockets 
are mainly attributed to d^z and dy^ orbitals respectively {g, 7|. Since the a band contains 
both odd and even components (named as and a„ respectively hereafter), the d^z with 
even symmetry alone cannot account for the observation. Moreover, the /3 band is largely 
suppressed in a geometry, and thus it must be mostly even, while theory predicted it to be 
odd (dyz). Similarly around M, the calculations suggest that the large and small electron 
pockets along T — M are made of the d^y and dyz orbitals respectively, both of which are odd. 
But in experiments, the large (7') and the small (7) pockets show even and odd symmetry 
respectively. 

There is another important constraint put by the experiment. Namely in the n geom- 
etry, the intensities of the same band on opposite sides of F show a strong asymmetry of 
R=I{kx+)/ I{kx-) ~ 5 in Figs. |2]^d) or|2]^e). To understand this. Fig. HJshows the calculated 
photoemission matrix elements in the vr geometry for the dxz , and dx2-y2 contributions to 
the photoemission matrix element at kx = iO.lsA"^, and kz in the range of to 5A~^. We 



have neglected the possib 



e complication from the kz dispersion, since it is not so prominent 



in the non-SDW state [ll,|l9i. To fit the experimental intensity asymmetry, the final state 
kz is about OAA^^ for the dxz case [Fig. Hl^a)]. The photoemission matrix element of d^^ is 
comparable with that of d^z, but has the opposite asymmetry at low kz- On the other hand. 
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FIG. 5: (color online) APRES data around T along the F-X direction, (a) The momentum cut 
for the presented data is marked with a thick solid line in the Brillouin zone. The symmetries of 
the five Sd-orbitals are illustrated, (b) Photoemission intensity and (c) its second derivative with 
respect to energy taken with vr geometry, (d) the EDO's for the data in panel b. (e) the MDC's 
and the fitting results for the data in panel b. (f, g, h and i) are the same as in panel b, c, d and 
e respectively, but taken with a geometry. 

the matrix element of d^^ can fit the experimental asymmetry at = 4AA~^ [Fig. lU^b)], 
where d^z contributes almost equally strong intensities to both sides. Thus, the asymmetry 
of our spectrum demands that either single d^i or d^z orbital should dominate the bands 
near F. While our current experimental result can not further pin down the exact orbital, 
considering the fast in-plane dispersion in the F — M direction and the known band structure 
calculations ^, ISj, we conclude that there is little contribution of dz^ to the Fermi surface 
near F. The photoemission cross-section of dx2_y2 is two order of magnitude smaller than 
the other two near F [Fig. ID^c)]. Therefore, although there may be strong mixing of d^2_y2 
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orbitals, it contributes negligibly to the measured photoemission intensity. This behavior of 
drr 2 y2 has been observed in cuprates, where very little spectral weight can be observed near 
r 20|. Because the intensity of the component is even somewhat stronger than that of 



the (5 band [Fig. E^e)], the /3 band should not be composed of single d^z or d^^^y^ orbitals, 
but rather strong hybridizations of both dxz and dx2_y2 orbitals, and the component 
contains more dxz state than the /3 band. On the other hand, our calculation shows that 
for low kx, the dxy orbital contributes little to photoemission intensity compared with the 
dyz orbital [FigsJU^d-e)], therefore, the component contains sizable dyz plus possible dxy 
orbital (Table. I). 

To further determine the orbital characters for the bands around F, we rotate the mirror 
plane by 45° with respect to z-axis, to the F-X direction. Photoemission intensity along the 
F-X direction in two geometries are shown in Fig. |5l The polarization dependence of the 
a and (3 bands resembles that in Fig. [2l The dxy and dz2 orbitals are even, and the dx2-y2 
orbital is odd in this setup. Although dxz and dyz orbitals are not symmetric as shown 
in Fig. |5]^a), their hybridized orbitals, '^'"^^"^ and '^'""'^"^ , have even and odd symmetry 
respectively. Therefore, along the F-X direction, the outer pocket should have dxy, dz2 or 
dxz+dyz Qj,]^j^g^j characters. The inner pocket should mainly be the mixture of both even 
and odd orbitals. Similar polarization dependence behaviors are observed when the mirror 
plane is rotated by another 45° to the yz plane, as expected from the 4-fold symmetry of 
the system. These indicate that the dxz portion in a.,^ and /3 along the F — M direction is 
evolving to '^"^'^y" in the T — X direction, and to dyz in the F — M direction; while the 
sizable dyz portion in evolves to '^"'^^^^ and to dxz respectively. 

We summarize our experimental orbital assignments of the bands in Table-I. A few re- 
marks are listed as follows: (1) Strong polarization dependence has been reported recently 
by several ARPES groups jsH]. These experiments were done on the parent compounds of 
iron-based superconductors. The analysis of orbital characters from these experiments was 
not straightforward due to the presence of lattice distortion and magnetic ordering. Our re- 

nn 

suits show significant differences compared to those in the SDW states [9|, lUj as pointed out 
above. (2) Our results suggest that there are distinct orbital characters for the most bands 
near Fermi surface. This property is very important for constructing microscopic models and 
investigating physical properties such as disorder effects. (3) Our comprehensive analysis 
does not fully rule out the possible major presence of dz2 at Fermi surfaces, which is possible 
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TABLE I: The main orbital characters of the bands near Ep along the T — M direction. 
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when strong electron-electron correlation is included 2l| . (4) The current five-band models 
cannot explain our experiments, which is not entirely unexpected. For example, the Fermi 
sphere in the T — Z direction, which was predicted by band calculations and attributed 
to in BaFe2As2 compounds, has never been observed so far. The discrepancy between 
the theoretical band structures and our experimental results suggests that electron-electron 
correlation may play an important role. 

To summarize, we have carried out a systematic investigation of orbital characters in 
BaCoo.i5Fei.85As2 and have analyzed all the possible orbital characters of each bands based 
on strong polarization-dependent photoemission response of the low-lying electronic struc- 
ture. We find that most bands near Fermi surface have distinct orbital characters. Our 
explicit conclusions cannot be understood in the current band calculations, and put strong 
constraints on the theory of iron-based superconductors. 
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NO.2006CB921300), STCSM of China, the NSF of US under grant No. PHY-0603759, and 
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X. Y. Cui and Prof. L. Patthey at SLS. 
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